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c 2 ^ ' The procedure of extracting the ground-state parameters from vacuum-to-vacuum and vacuum- 

^ to-hadron correlators within the method of sum rules is considered. The emphasis is laid on the 

. crucial ingredient of this method — the effective continuum threshold. A new algorithm to fix this 

(N ■ quantity is proposed and tested. First, a quantum-mechanical potential model which provides the 

$H only possibility to probe the reliability and the actual accuracy of the sum-rule method is used as 

a study case. In this model, our algorithm is shown to lead to a remarkable improvement of the 
' accuracy of the extracted ground-state parameters compared to the standard procedures adopted 

I in the method and used in all previous applications of dispersive sum rules in QCD. As a next 

, step, it is demonstrated that the procedures of extracting the ground-state decay constant in the 

potential model and in QCD are quantitatively very close to each other. Therefore, the application 
' 1^ ' ' of the proposed algorithm in QCD promises a considerable increase of the accuracy of the extracted 

^^P^- hadron parameters. 

^ . 1. INTRODUCTION 

The method of dispersive sum rules [l|, Q is one of the widely used methods for obtaining properties of the ground- 
" ' I state hadrons in QCD. The method involves two steps: (i) one calculates the relevant correlator in QCD at relatively 

: small values of the Euclidean time; (ii) one applies numerical procedures prompted by quark-hadron duality in order 
\^ ' to extract the ground-state contribution from this correlator. These numerical procedures cannot determine a single 

. value of the ground-state parameter but should provide the band of values containing the true hadron parameter with 
y—i ' a flat probability distribution within this band flj. The width of this band is a systematic, or intrinsic, uncertainty 

, of the method of sum rules. 

' It is often claimed that the standard procedures of the method of sum rules allow one to obtain hadron parameters 
I with about 10% accuracy for normal and exotic hadrons. In many cases this is true. However, there are some 
T-H . examples where the difference between the sum-rule predictions and the experimental results or between the results 
^ ' from different versions of sum rules turn out to be quite unexpectedly much larger than 10%: as the first example, recall 
• I— I , an evident conflict between the prediction for the D*Dit coupling from light-cone sum rules Q and the subsequent 
' measurement of CLEO as the second example, mention the incompatibility of the predictions for the pion form 
^ , factor at intermediate momentum transfers from light-cone sum rules [F] , local-duality sum rules |6| , and sum rules 
with non-local condensates Q- A posteriori, these issues may be clarified. The real problem is, however, that the 
large errors — well exceeding a comfortable 10% level — could not be foreseen on the basis of the standard criteria 
adopted in the method of sum rules in QCD. 

An unbiased judgement of the reliability of the extraction procedures adopted in the method of sum rules may be 
acquired by applying these procedures to problems where the ground-state parameters may be found independently 
and exactly as soon as the parameters of the theory are fixed. Presently, only quantum-mechanical potential models 
provide such a possibility. 

A simple harmonic-oscillator (HO) potential model Q possesses two essential features of QCD — confinement and 
asymptotic freedom Q — and has the following valuable features: (i) the bound-state parameters (masses, wave 
functions, form factors) are known precisely; (ii) direct analogues of the QCD correlators may be calculated exactly. 
(For a discussion of many aspects of sum rules in quantum mechanics we refer to p^l - flsj V 

Making use of this model, we have already studied the extraction of ground-state parameters from different types 
of correlators: namely, the ground-state decay constant from two-point vacuum-to- vacuum correlator the form 
factor from three-point vacuum-to- vacuum correlator [isl] , and the form factor from vacuum-to-hadron correlator JlG^] . 
We have demonstrated that the standard procedures adopted in the method of sum rules do not work properly: for all 
types of correlators the true value of the bound-state parameter was shown to lie outside the band obtained according 
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to the standard criteria. These results gave us a sohd ground to claim that also in QCD the actual accuracy of the 
method turns out to be worse than expected on the basis of applying the standard criteria. 

We have realized that the main origin of this difficulty of the method lies in an over-simplified model for the 
hadron continuum, which is described as a perturbative contribution above a constant Borel-parameter-independent 
effective continuum threshold. We have introduced the notion of the exact effective continuum threshold, which 
corresponds to the true bound-state parameters: in the HO model the true hadron parameters — decay constant and 
form factor — are known and the exact effective continuum thresholds for different correlators may be calculated. 
We have demonstrated that the exact effective continuum threshold (i) is not a universal quantity and depends on 
the correlator considered (i.e., it is in general different for two-point and three-point correlators), (ii) depends on the 
Borel parameter and, for the form-factor case, also on the momentum transfer. 

Moreover, we have shown that the "Borel stability criterion" combined with the assumption of a Borel-parameter- 
independent effective continuum threshold leads to the extraction of the wrong bound-state parameters. 

In recent publications [l7l - [l9j we formulated a new algorithm for extracting the parameters of the ground state. 
The idea proposed in these papers is to relax the standard assumption of a Borel-parameter-independent effective 
continuum threshold and allow for a Borel-parameter-dependent quantity. 

In the present paper we develop this idea and provide details of its application to the extraction of hadron parameters 
from dispersive sum rules. We show that (i) the new algorithm leads to improvements — and for the case of the 
form factors to significant improvements — of the actual accuracy of the extracted bound-state parameters from all 
kinds of correlators in the potential model; (ii) the procedures of extracting the ground-state parameters from the 
truncated operator product expansion (OPE) for the two-point function in potential models and in QCD exhibit not 
only qualitative but also quantitative similarities. 

The paper is organized as follows: Section [2] presents the application of our new algorithm to the extraction of the 
ground-state parameters from correlators in a quantum-mechanical potential model. Quantum-mechanical analogues 
of all correlators used in QCD for the extraction of ground-state parameters (decay constants and form factors) are 
discussed. Section [3] compares the extraction of the ground-state decay constant in QCD and in quantum-mechanical 
potential models. Section 2] presents our concluding remarks. 

2. HARMONIC-OSCILLATOR POTENTIAL MODEL 

We consider a nonrelativistic Hamiltonian H with a HO interaction potential V(r), r = |r|: 

H = Ho + V(r), Ho ^ pV2"^, V{r) = (2.1) 

In this model, all characteristics of the bound states are easily calculable. For instance, for the ground state (g), one 
finds 

3 / 777 / . ; \ 3/2 

E,^-u;, i?g^|^g(r = 0)p = (--) , i^g(g)=exp(-<jV4ma.), (2.2) 
where the elastic form factor of the ground state is defined in terms of the ground-state wave function ijjg as 

FM = (V'gk(q)IV'g) = J rf3fcV4(k)Vg(k-q) - J dV|^g(r)pe*'^^ |q|, (2.3) 

and the current operator J(q) is given by the kernel 

(r'l J(q)li-) - exp(iqr)5(3)(r _ r'). (2.4) 

A. Polarization operator 

The quantum-mechanical analogue of the Borelized polarization operator has the form 

H(r) = (r/-0|cxp(-i/T)|r, = 0). (2.5) 
For the HO potential, the analytic expression for H(r) is well-known [§]: 



n(r) 



( — ) . ^,.3/2 - (2-6) 

[2 smh(cL>J )J ' 
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Apart from the overall factor, n(T) is a function of one parameter loT. 
Let us define the average energy of the polarization function as 

En{T) = -aTlogn(r) = -c^coth(wr), dr = (2.7) 

At T = both n(T) and Eii{T) diverge. 

Making use of the "hadronic" bound states of the model, we obtain the following representation for n(T): 

nhadr(T) = ng(T) + ncxcitcd(r), ng(T) = i?g exp {-e^t) . (2.8) 

For large values of T the contributions of the excited states to the correlator vanish and therefore Eu{T) — ^ Eg for 
T — > cx). The deviation of the energy from Eg at finite values of T measures the "contamination" of the correlator by 
the excited states. 

The OPE series is the expansion of n(T) at small Euclidean time T: 

3/2 1 o„o 19 



For n(T) this expansion is equivalent to the expansion in powers of the interaction tu. The first term in (j2.9l) . no(r), 
does not depend on the interaction and describes the free propagation of the constituents. The rest of the series 
represents power corrections: 

npower(r) = noPE(r) - no(r). (2.10) 

no(T) may be written as the spectral integral 

no(r) = dze-'^po{z), po{z) = ^^V^: (2-11) 
with po {z) the spectral density of the one- loop two-point Feynman diagram of the nonrelativistic field theory [T^ | . 

B. Vertex function 

The basic quantity for the extraction of the form factor in the method of dispersive sum rules is the correlator of 
three currents [2] . The analogue of this quantity in quantum mechanics reads p3l [l5| 

T{E2,Ei,q) = {vf=0\{H-E2)-'J{q){H-Ei)''\r,^0), g=|q|, (2.12) 

[with the operator J(q) defined in (|2.4p ] and its double Borel (Laplace) transform under Ei ti and E2 -> T2 

r(T2,ri,g) = (r^ =0|G(r2)J(q)G(Ti)|r, -0), G(r) = exp(-iJT). (2.13) 

For large ti and T2 the correlator is dominated by the ground state: 

r(T2, n, g) ^ |V'g(r = 0)|2e-^-(^i+^^'Fg(g2) + . . . . (2.14) 

Let us notice the Ward identity which relates the vertex function at zero momentum to the polarization operator: 

r(T2,ri,g = 0) =n(Ti +T2). (2.15) 

This expression follows directly from the current-conservation relation 

J(q = 0) = l. (2.16) 

Let us consider the vertex function for equal times ti — T2 — ^T, which in the HO model explicitly reads 

r{T,q) = n(r)exp L-^tanh f'^)) . (2.17) 



Amui \ 2 J 

The correlator F is a function of two dimensionless variables luT and , q = q/^/rmv 
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The corresponding average energy is defined as follows: 

EriT,,) . -5.1ogr(r,,) . ^.coth(.T) + £^j^^^;^^y (2.18) 

The contribution of the ground state to the correlator, Fg, in the HO model has the form 

r^{T,q)^R,exp{-E^T)F,{q'). (2.19) 

The relative contribution of the ground state to E-r{T, q) and r(r, q) are shown in Fig. [TJ independently of the value 
of g, E{T, q) — >■ i?g and F(T, q) — ^giT, q) for large T. Notice, however, that the saturation of the correlator by the 
ground state is delayed to later T with the growth of q. 

Let us now construct for F(r, q) the analogue of the OPE as used in the method of three-point sum rules in QCD. 
First, we expand F in powers of and obtain 

oo 

F(T,(?) = ^F2„(r,g)c.2„_ (2.20) 

n=0 

This expansion is the analogue of the expansion of the vertex function in terms of the nonlocal condensates in QCD 
[2l| . Explicitly, for the lowest terms one has 



ro(r,,) . 1^ 



.... 

The first term, Fq, corresponds to free propagation and does not depend on the confining potential. It may be written 
as the double spectral representation [22| 

/•oo 

Fo(r,g) = / dzidz2e-^''^e'^''^Ao{zi,Z2,q), (2.22) 
Jo 

Mz,,Z2,q) = ^^((-i+-2-|^) -4z,z,<oj. 

The Ward identity in nonrelativistic field theory relates to each other the three-point function at zero momentum 
transfer and the two-point function and leads to 

lim Ao(z, z', q) = S{z - z')po{z). (2.23) 

Power corrections are defined as the difference between the exact expression and the free diagram: 

roPE(r, q) = Fo(T, q) + Fpowcr(r, q). (2.24) 

In order to obtain the analogue of the OPE for F in terms of local condensates, one should further expand Fpowcr(2^, <z) 
(i.e., F2„(r, g), n > 1) in powers of T. 
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C. Vacuum-to-hadron correlator 



In order to apply the method of sum rules to hadron form factors at intermediate and larg e momentum transfers 
a Borelized vacuum-to-hadron amplitude of the T-product of two quark currents is used |2C| . The analogue of this 
quantity in quantum mechanics has the form jT6| 

A{T, q) = {r = 0|G(r) J(q)IV'g), G{T) = exp(-iJT). (2.25) 

The spectral representation for A(T, q) may be written in the form 



A{T,q) 

where Q 



j d^kG{k^,T)iP^{\i~q)^ dz ejipi-zT)pAiz,T,q), z = k^/2m, (2.26) 



G{k',T) ^ (r = 0|G(T)|k> - \ wV,i3/2 ^^P f-;^^ tanh(c.r)) . (2.27) 

(27r)-^/^ [cosh(a;T)J'^'"^ \ 2mw / 

Explicit expressions for pa{z, T, q) and A{T, q) may be obtained using i/jg of Eq. (|2.2p : for instance, for A{T, q) one finds 

AiT, q) = /R^exp i^E^T) exp (^-^ (l - g-^-^)^ . (2.28) 

The function A{T, q) depends on two dimensionless variables (f' ~ q^ jmu) and ujT . 

Expanding in Eq. (I2.26P the Green function G{k^,T) in powers of uj generates the analogue of the QCD twist 
expansion for the amplitude A{T,q) (see [l^ for details). 

The ground-state contribution to the correlator has the form 

A^{T,q) = ^,cxp{~E,T)F^iq), (2.29) 

so one finds 

^MZ4^expf-^e--^V (2.30) 

A{T, q) ^ \ Amuj ) ^ ' 

Notice the following features of the correlator A{T,q) (see Fig. [5]): 

(i) At large T, the ground state provides the dominant contribution to A, similar to the case of vacuum-to-vacuum 
correlators: 

A{T,q) ^ Ag_{T,q) for T ^ oo. (2.31) 

(ii) At small q, because of current conservation, the ground state dominates the correlator for all T. This is a specific 
feature of A which arises due to the choice of the initial state. Thus, compared with vacuum-to- vacuum correlators, 
the correlator A "maximizes" the ground-state contribution. 
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D. Dual correlators and effective continuum thresholds 

The sum rule for each of the correlators described above is an expression of the fact that the correlator calculated in 
terms of the constituent degrees of freedom (quarks and gluons) and the bound-state degrees of freedom (hadrons) are 
equal to each other. In order to isolate the ground-state contributions from the sums over hadron spectra, one invokes 
quark-hadron duality which assumes that the excited states are dual to the high-energy regions — i.e., the regions 
above some effective continuum thresholds — of the spectral representations for the correlators. The correlators to 
which the corresponding cuts have been applied are referred to as the dual correlators. The relations between the 
ground-state contributions and the dual correlators take the following form: 

i?ge-^-^ = nduai(T,05f(r))= J dze-'^poiz)+Il„iT), (2.32) 



Fg{q)Rge'''-^ = rduai(r, g, z^^ (T, g)) = J dz^ J e~^("i+"^)^Ao(zi, ^2, g) + rpower(r, g), (2.33) 



Fg(g)/R;e-^-^ = Aduai(T, q, 2,^^ (T, q)) = / dz e-'^ pA{z,T,q). (2.34) 

As consequence of (|2.15p . npowGr(T) — rpowcr(I', 9 = 0) and 

^cff(r) = z[ff(r,q = 0). (2.35) 

Then, due to the Ward identity (|2.23p . the relations (|2.32p and (|2.33p yield the correct normalization of the form 
factor Fg^{q = Q) = 1, as required by current conservation. 

Let us take advantage of knowing the exact ground-state parameters in the potential model and calculate the exact 
effective continuum thresholds for all the correlators. These exact thresholds are obtained by solving the equations 
above for the exact ground-state parameters on the l.h.s.; the exact thresholds make the quark-hadron duality relations 
(|2.32p - (j2.34p exact. Clearly, z^^actC^i? = 0) = z^^^^^^{T). Figure |3] depicts the exact effective continuum thresholds 
for r(T, q) and A{T,q). As is obvious from Fig. [31 the effective continuum thresholds do depend considerably on 
both T and q. Moreover, for the same theory (i.e., the same potential and the same values of its parameters) the 
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exact effective thresholds for r(T, q) and A(T, q) are rather different from each other and have very different T- and 
q-dependencies. 

It might be useful to notice that the dependence of the effective thresholds on the Borel parameter, being unusual 
for spectral representations, does not contradict any properties of field theory: the dual correlators are hand-made 
objects and they cannot be written as, e.g., the T-product of local currents; thus the dual correlators have completely 
different analytic properties than the "usual" correlators in field theory.^ 

To summarize, the exact effective continuum thresholds (i) depend on T and on q; (ii) are not "universal", i.e., for 
one and the same ground state, the effective thresholds in different correlators strongly differ from each other. As a 
consequence, making use of the effective continuum threshold obtained, e.g., for Ilduai as effective thresholds for the 
dual correlators Fduai or A^^ai, as often done in practical applications of the sum-rule method in QCD, has no ground 
and may lead to very inaccurate estimates for the form factors. 

On the practical side, a good or a bad extraction of the ground-state parameters depends on our capacity to find 
a reasonable approximation to the exact effective continuum threshold. 



E. Numerical analysis 



Let us consider a restricted problem where the energy Eg of the ground state and its wave function at the origin, 
y/R^, are known, and try to extract its elastic form factor from the sum rules p.33p and ()2.34p . 

A. First, according to [l| we should determine the working interval of T — the "Borel window" — where the method 
may be applied to the extraction of the ground-state parameter: 

(i) The upper boundary of the T-window is obtained from the condition that the truncated expansion gives a good 
approximation to the exact correlator. Since in the HO model the correlators are known exactly, the upper boundary 
is T = oo. However, to be close to realistic situations, where only a limited number of higher-order or higher-twist 
corrections is available, following our results reported in [isl . [l6j we define the upper boundary of the window as 
wT < 1.2 for T{T,q) and ujT < 0.8 for A{T,q). 

(ii) The lower boundary of this T-window is determined by the condition that the ground state gives a "sizable" 
contribution to the correlator. If we require the ground-state contribution to exceed, say, 50%, the window for r(r, q) 
closes already for g ~ 2 (Fig. [T})), whereas for A{T, q) it closes for g ~ 4 (Fig. [2}3). 

We shall see, however, that our algorithm will allow us to extract the form factor, although with a worse accuracy, 
also in the region where the ground-state contribution to the correlator remains well below 50%. Thus, our algorithm 
opens the possibility to study also the region of larger momentum transfers. 

We choose the "windows" as follows: 

0-7<a;r<1.2 for T{T,q), (2.36) 
0.2<a;T<0.8 ioi A{T,q). (2.37) 

B. Second, we must choose a criterion to fix the effective continuum threshold Zcff (r, q). We proceed in the following 
way: We consider a set of T-dependent Ansiitze for the effective continuum threshold 

n 

^iff (r,g)=E#^(9)(^T)^'. (2.38) 

(The standard procedure adopted in all sum-rule applications in QCD is to assume a T-independent quantity. Obvi- 
ously, this option is also included in our analysis.) Now, at each value of q we fix the parameters on the r.h.s of (|2.38p 
as follows: we calculate the dual energy 

ELAT,q) - -^logFduai(T,<?,zfff(T,q)), (2.39) 
E^uAT,q) = -^logAduai(T,q,z,1j(T,g)), (2.40) 



The dual correlator has a relatively simple form in the Borel-parameter space; the dual correlator in the momentum space is a nonlocal 
object. This is not strange: the dual correlator should reproduce a single ground-state pole by the dispersion integral over a finite energy 
region. Therefore, the dependence of the effective threshold on T is its inherent property required by the analytic properties of the dual 
correlator. 
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(a) (b) 

Fig. 4: The dual form factor -Fduai(g) extracted from (a) the correlator r{T,q) and (b) the correlator A{T,q). Red symbols 
(n = 0) correspond to the standard T-independent Ansatz for ZcS, green (n = 1) to the linear Ansatz, blue {n = 2) to the 
quadratic Ansatz, and black (n — 4) to the quartic Ansatz. 



for the T-dependent of Eq. (|2.38p . Then we evaluate i^duai {T, q) at several values oiT ~ Ti [i — 1, . . . , N , where 
can be taken arbitrarily large) chosen uniformly in the window. Finally, we minimize the squared difference between 
i^duai and the exact value E^: 



N 

,2 _ 



X 



^EK;P(^-«)-^J'- (2-41) 



Recall that the exact effective continuum thresholds for the two correlators V{T,q) and A{T,q) have rather different 
T-depcndcncies, see Fig. [31 Therefore, one should expect that the expressions (|2.38p for the corresponding effective 
thresholds will be different. 

The results for the form factor obtained from r(T, q) and from A{T^ q) upon optimizing the parameters of z^g and z^g 
according to (|2.4ip are shown in Fig. U) 



1. Form factor from T{T, q) 



Let us first consider the extraction of the form factor from the correlator T{T,q). 

Within the standard T-independent approximation, one extracts the form factor in the region q < f .5 with better 
than 10% accuracy, while the accuracy decreases rather fast for higher q (see Fig. |4^). The dramatic increase of the 
error at g > 2 is related to the fact that the contribution of the ground state to the correlator decreases rapidly with 
q in the given T-window 0.7 < ujT < 1.2 (see Fig.[T|). 

The real problem is, however, that the magnitude of the error cannot be obtained on the basis of the standard 
criteria adopted in the method: e.g., at g = 2, the variation of the extracted form factor in the window is only about 
2% (see Fig. Eb) which mimics an accurate extraction of the form factor, whereas the actual error comprises 25%. 
Thus, within a T-indepcndent Ansatz for the effective threshold, the error of the form factor cannot be determined 
from the variation of the dual form factor in the window. 




Fig. 5: The correlator T{T,q): the fitted dual energy EduaiiT,q) and the corresponding dual form factor Fduai{T,q) for q — 2. 
Red lines (n = 0) correspond to the standard T-independent Ansatz for Zcfi, green lines (n = 1) to the linear Ansatz, and blue 
lines (n — 2) to the quadratic Ansatz. 
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Allowing for the T-dependent effective threshold leads to a considerable improvement of the accuracy: both for the 
linear and the quadratic Ansatze for the threshold, the dual energy reproduces the known energy rather accurately 
in the window, and the dual form factor reproduces the exact form factor rather well. Noteworthy, the actual error 
of the extracted form factor may be obtained as the spread between the dual form factors obtained for the linear and 
the quadratic Ansatze for z^g. 



2. Form factor from A{T, q) 



Let us now consider the extraction of the form factor from the correlator A{T,q). 

The standard T-independent Ansatz works well below q < 1.5 — 2, but for larger q the accuracy of the extracted 
form factor decreases very fast, see Fig. |4l The main problem here is the same as for the case of T{T,q): the large 
error in the extracted form factor cannot be guessed on the basis of the standard Borel stability criterion: as seen in 
Fig. [HI the variation of the form factor at g = 2 in the window comprises only 2%, whereas the actual error exceeds 
10%. 

Allowing for T-dependent approximations for leads to considerable improvements: as shown in Fig. |4l the form 
factor may be extracted with a reasonable accuracy (better than 15%) up to g = 5 with the quadratic and quartic 



Ansatze, z^{T, q) and z^g{T, q) (the results obtained for n ~ 1 are close to the results for n = 0, and the results for 
n = 3 are close to those for n = 2). 

This high accuracy is achieved in spite of the fact that (i) the relative contribution of the ground state to the 
correlator falls down considerably below 50% in the T-window, and (ii) the form factor at q — 5 falls down compared 
to its value at g = by almost three orders of magnitude. 

Notice again the typical feature of the extraction procedure: both for n = 2 and rt = 4 both the dual energy 
Edua\{T,q) and the dual form factor ^duai(r, g) are extremely stable (better than 0.1%) in the window for all q < 5 
(Fig.[S]shows the results for q = 2). In spite of this stability, the actual error of the extracted form factor for n = 2 and 
n = 4, separately, is found to be much larger: e.g., at q = 5 this error comprises 10-15%. Again, this error could not 
be guessed on the basis of the standard Borel stability criterion, which thus does not work also for vacuum-to-hadron 
correlators. 

Nevertheless, we have seen clear improvements in the results of the extraction procedure as soon as one goes beyond 
the assumption of a T-independent effective continuum threshold: First, the actual accuracy turns out to be (much) 
better. Second, the T-dependent Ansatz allows one to extract the form factor in a broader range of the momentum 
transfer. Finally, the spread between the results obtained with the quadratic and the quartic Ansatze for Zcs provides 
the actual error of the extracted form factor. 



.(4) 
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Fig. 6: The correlator A{T,q): the fitted dual energy Edua.\{T,q) and the corresponding dual form factor Fduai{T,q) for q — 2. 
Red lines (n = 0) correspond to the standard T-independent Ansatz for ZcH, blue lines (n = 2) to the quadratic Ansatz; and 
black lines (n — 4) to the quartic Ansatz. 



3. POTENTIAL MODEL VS. QCD 



In the previous section we have shown that the concept of the T-dependent effective continuum threshold and the 
way of fixing this quantity by reproducing the bound-state energy has led to a dramatic increase of the actual accuracy 
of the extracted form factor in the potential model. 

Does this new approach to the extraction of the properties of ground-state hadrons promise equivalent improvements 
in QCD? We give a positive answer to this question [23|. Figure [7] compares the extraction of the ground-state decay 
constant from the two-point function in the potential model (left column) and in QCD (right column). For the 
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Fig. 7: Dual mass and extracted dual decay constant vs. Borel parameter in the potential model (left column) and in QCD 
(right column). In QCD, the Borel parameter r has dimension GeV~^ . 



potential model, we consider the interaction consisting of a confining HO part and a Coulomb part and make use of 
the OPE which includes the radiative corrections up to order and several lowest-order power corrections. This 
corresponds to the best known OPE for the two-point function of heavy-light quark currents [l^j which we use for 
QCD. Details of this analysis including the specific values of the numerical parameters can be found in 23]; they are 
not essential for our purpose here. 

From Fig. [7] it becomes evident that — with respect to the extraction of the ground-state parameters — there are 
no essential differences between QCD and quantum mechanics: as soon as the parameters of the Lagrangian are fixed, 
and the truncated OPE is calculated with a reasonable accuracy (taking into account also the relevant choice of the 
renormalization scale in QCD), the extraction procedures are very similar. 

At first glance, this similarity might look surprising since the structure of bound states in quantum-mechanical 
potential models and in QCD are rather different. However, the method of dispersive sum rules does not make use 
of the details of the ground-state structure. What really matters for extracting the ground-state parameters by this 
method is (i) the way one implements the quark-hadron duality and (ii) the structure of the OPE for a given correlator. 
Therefore, it should not be surprising that the extraction procedures in potential models and in QCD are similar, as 
soon as quark-hadron duality is implemented in both theories in the same way by a cut on the spectral representation 
for the correlators. 

In the potential model, the actual value of the decay constant belongs to the band provided by the linear and 
the quadratic Ansatze for the effective continuum threshold. It seems very likely that the same happens in QCD, 
at least this may be taken as a realistic conjecture. Perhaps even more important is that studying the r-dependent 
parameterizations for the effective threshold allows one to probe the magnitude of the intrinsic, i.e., systematic, 
uncertainty of the extracted decay constant. Interestingly, the systematic uncertainties of the ground-state parameters 
are probed solely by the OPE, without making use of the information about the excited states of the hadron spectrum. 
The intrinsic uncertainty should be taken into account along with the "usual" uncertainty obtained in the method of 
QCD sum rules by the variation of the input QCD parameters; while the latter is going to decrease when increasing the 
accuracy of these parameters, e.g., mhijnb), the systematic uncertainty does not depend on the input QCD parameters 
and remains at the same level. 



4. CONCLUSIONS 



We studied the extraction of the ground-state form factor from the vacuum-to-vacuum and vacuum-to-hadron 
correlators. Let us summarize the main messages of our analysis; 

• The knowledge of the correlator in a limited range of relatively small Euclidean times T (i.e., large Borel masses) 
is not sufficient for the determination of the ground-state parameters. In addition to the OPE for the relevant 
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correlator, one needs an independent criterion for fixing the effective continuum tfiresliold. 

• Assuming a T- independent (i.e., a Borel-parameter-independent) effective continuum tfiresfiold, tlie error of tlie 
extracted dual form factor JduaK?", q) turns out to be typically much larger than the variation of F(T, q) in the 
Borel window. Therefore the standard procedures based on the Borel stability adopted in the method of sum 
rules do not allow one to obtain realistic error estimates for the extracted form factor; eventually the actual 
error of the extracted form factor may happen to be much larger. 

• In those cases where the ground-state energy (mass) is known, considerable improvements may be achieved by 
allowing for a T-dependent effective continuum threshold and finding its parameters by minimizing the deviation 
of the energy of the dual correlator from the known value of the ground-state energy, Eq. (|2.41|) . First, the actual 
accuracy of the extracted form factor is considerably improved. Second, the spread between form-factor values 
obtained with different T-dependent Ansatze for Zeff provides a realistic error of the form factor. 

• We have demonstrated the striking similarity — not only qualitatively but, more importantly, also quantitatively 
— between the extraction procedures in the quantum-mechanical potential model and in QCD. This similarity 
is the consequence of the same way of implementing the quark-hadron duality in both theories. The established 
similarity strongly suggests that the results on the extraction of ground-state form factors obtained in potential 
models have direct implications for the corresponding analyses in QCD and should be taken very seriously. 

Although a rigorous — in the mathematical sense — control over the systematic errors of the extracted hadron 
parameters remains unfeasible, the proposed modifications related to the Borel-parameter-dependent effective 
continuum threshold promise a considerable increase of the actual accuracy of the method of sum rules and 
opens the possibility to obtain realistic estimates for the systematic errors. 
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